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I. INTRODUCTION TO THE INITIAL-BOUNDARY VALUE PROBLEM IN GENERAL RELATIVITY 

This article has been written to celebrate Mario Castagnino's seventy fifth birthday. Knowing his particular interest 
in fundamental problems which has marked his scientific prolific career we are confident he will share our curiosity 
in this subject. We consider a manifold M of the form M = [0, T] x £, where [0,T] is a time interval and S a 
three-dimensional manifold with smooth boundaries <9£, see Fig. [1] We want to answer the following question. What 
data must be given at the initial surface £ := {0} x £ and the boundary surface T := [0,T] x d£ such that there 
exists a unique solution (M, g) of Einstein's vacuum field equations which depends continuously on the data? From 
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FIG. 1: The geometrical setup for the initial-boundary value formulation. 
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the work of Mme Choquet-Bruhat[l| in the fifties we know that at So wc must give a pair of symmetric tensor fields 
(h a b,k a b), the first being positive definite (and so a metric), which obeys the Einstein constraint equations. They 
are called the initial data sets for Einstein's equations. They give a unique solution to Einstein's equations on the 
domain of dependence -D(Eo) of £0 (which depends on the data given), namely a spacetimc (D(£o),g), such that 
ip*g = h,tp*C n g = k, where n denotes the unit normal to £0 m D(T,q) and tp : £0 ~~ ^ D(T,q) the inclusion map. This 
solution depends continuously on its initial data. This setting is called the initial value problem. Furthermore, if two 
initial data sets are related by a diffeomorphism in £0, then the unique solutions they produce are also related by a 
diffeomorphism, this time on the maximal development [H of £o- We call this property the geometric uniqueness (or 
at least a version of it) of the initial value problem. 

To show the above assertion is not an easy nor a direct task since Einstein's equations are geometrical and so they 
are subject to the diffeomorphism freedom. First, the field equations need to be recast into a system of evolution 
equations with constraints. Second, a gauge for which the first set of equations is strongly hyperbolic needs to be 
found. Third, the theory of strongly hyperbolic equations is used in order to show that a unique solution in this 
gauge class exists and depends continuously on the initial data. Finally, geometric uniqueness must be shown. This 
program was first accomplished using harmonic coordinates. 0, Q 
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Contrary to the initial value problem the data at the boundary T is not unique. Here, we are thinking about the 
case where T represents a time- like surface (with respect to the unknown metric g). Therefore, what can be specified 
at the boundary is, essentially, the incoming radiation which could be given as an absolute quantity, zero for all times, 
say, or which could be given as a function of the outgoing radiation, as is the case for a mirror in electromagnetism, 
for instance. A further complication which makes the problem difficult to tackle arises from the absence of a reliable 
local definition for gravitational radiation in general relativity. Ultimately, this problem is related to the fact that 
the initial-boundary value problem (IBVP) we are considering appears as a means for having an artificial boundary 
for numerical calculations, and not an IBVP with physical boundaries which would, presumably, dictate the correct 
boundary conditions. This problem is, in fact, connected to two other ones: the first one is that in most formulations 
some of the constraint variables, that is, the quantities which yield the constraint equations when set to zero, have 
nonzero propagation speeds with respect to the boundary. In particular, there are constraints propagating into the 
computation domain, and in this case the boundary data must be chosen with care to ensure that no constraint- 
violating modes enter the domain. The second problem, which conspires with the first one in making the task 
difficult, is that in most formulations the problem is characteristic. We shall provide a precise definition shortly, but 
roughly speaking it means that there are perturbations having zero speeds with respect to the boundary, that is, 
modes which propagate along the boundary. It is hard to control such perturbations at the boundary, and so it is 
difficult to formulate a well posed problem if such zero speed perturbations appear in the boundary conditions fixing 
the constraints. 

The IBVP for Einstein's vacuum field equations was solved by Friedrich and Nagy in 19990 based on a frame 
formalism with very precise gauge conditions tailored to the problem. Furthermore, their formulation also modifies 
the evolution equations near the boundary in such a way that the constraint variables propagate along the boundary. 
This feature implies that no constraint conditions arc needed at the boundary and one gets away with boundary 
conditions which only control the physically relevant field components. After this impressive work most attempts to 
solve the problem for metric based formulations failed, the reason being that several of the issues mentioned above 
still needed to be understood. 

The first breakthrough came in 2006 when Kreiss and WinicourQ realized two very important points. The first 
point is that the second-order wave equations, when viewed as a certain pseudodiffcrential first-order operator, is 
not characteristic. This could be appliedJBJ to Einstein's field equations in harmonic coordinates, in which case 
the evolution equations reduce to a quasilinear system of wave equations for the metric components. Second, they 
realized that it was possible to specify constraint-preserving boundary conditions for this system similar to Sommerfeld 
boundary conditions, i.e. non-incoming radiation conditions, but with a certain coupling obeying a hierarchy which 
allowed to prove strong stability (defined below). More recently, the above results were also established using the 
usual energy estimates based on integration by parts. @ Furthermore, it was possible to conclude that the boundary 
conditions imposed were actually of the maximal dissipative type0, Q for a specific, non-standard class of first- 
order symmetric hyperbolic reductions for which the system is not characteristic. This also opens the interesting 
possibility to extend the result to discretizations of the equations using finite differences satisfying "summation by 
parts". [Ej 

The plan for the remainder of this paper is as follows. In the next section we discuss boundary conditions for 
strongly hyperbolic systems, introduce the concept of strongly stable problems and assert the main theorem on the 
subject together with two applications. In section IIIII we explore the boundary conditions for the linearized theory 
and analyze the geometric uniqueness problem. Finally, conclusions and open issues are discussed in section HVl 

II. STRONGLY STABLE IBVP 

Consider a first-order quasilinear system of the form 

d t u a = A ac {u, x, t)V c u p + B a (u, x, t), (1) 

where u a = u a (x,t), a = 1,2, ...N, are the unknowns, V c is a connection on M and A ac p(u, x, t) and B a (u,x,t) 
depend smoothly on their arguments. We say the above system is strongly hyperbolic if there exists a symmctrizcr, 
that is, a symmetric, positive definite matrix H a p = H(u, x, t, n) a p, depending smoothly on its arguments, such that 
H al A lc pn c is also symmetric for all (u, x, t) and all one-forms n c . The system is called symmetric hyperbolic if, in 
addition, the symmctrizer can be chosen independent of n. A quasilinear first-order system has a well posed initial- 
value problem which is stable under zeroth order perturbations (changes in B a (u,x,t)) if and only if it is strongly 
hyperbolic. Most physical problems describing wave propagation, i.e. electromagnetism, fluids, elasticity, general 
relativity, etc, can be cast into first-order symmetric or strongly hyperbolic systems, see for instance Ref. [Ill ] and 
references therein; in particular see Refs. [l2[ and [l3[ for covariant definitions. 
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Modulo the smoothness assumptions strong hyperbolicity is equivalent to asking that for each n c the principal 
symbol A ac p(u,x,t)n c has a complete set of eigenvectors with real eigenvalues. This implies that, locally, every high 
frequency solution is a combination of plane waves where the amplitude, frequency and wave vector of the wave is 
given by, respectively, the eigenvector, the eigenvalue and the one-form n c of the principal symbol. With the help 
of the symmetrizer it is possible to construct norms giving a priori energy estimates, a basic step in showing well 
posedness for the initial value problem. 

When boundaries are present, a systematic method 1 for constructing boundary conditions consists in analyzing the 
plane waves with normal incidence to the boundary, that is, in diagonalizing the boundary matrix 

N a (u,x,t) := A ac p(u,x,t)n c , 

where n c is an outgoing normal one-form to the boundary surface. 2 The field perturbations Su a can be written as 
a linear combination of the eigenvectors. The coefficients corresponding to positive, negative and zero eigenvalues 
are called the incoming, outgoing and zero speed fields, respectively. One must specify a condition for each incoming 
field. In contrast to this, the outgoing fields are determined by the evolution of the solution up to the given time. 
The zero speed fields are also predetermined by the evolution, but unfortunately in a weaker sense in comparison 
with the outgoing ones, and this is the main problem of the theory. The presence of these zero speed fields, that 
is the existence of zero eigenvalues of the matrix N, is usually referred as a problem with characteristic boundary. 
The boundary condition specifies data to the incoming fields, or to an appropriate linear combination of the in- and 
outgoing fields, as is the case for reflecting boundaries. However, the linear combination cannot include zero speed 
fields. Unfortunately, in theories with constraints, the boundary conditions required to enforce the constraints do 
usually contain such zero speed fields. 



A. Example: Maxwell's equations in the potential version 



In order to illustrate this problem we consider Maxwell's equations in the Lorentz gauge, in which case on obtains 
a set of wave equations for the components of the vector potential Ab- These equations are only coupled among them 
by the constraint the gauge imposes, namely the requirement for Ab to be divergence-free. For simplicity we just 
consider the theory in Minkowski spacetime, in which case 

V a V a A b = 0, V b A b = 0. 

Introducing the first-order derivatives of Ab as new independent fields, D a b := \? a Ab, the set of wave equations is 
reducible to the first-order system 

d t A b = D t b, d t D t b = d 3 D jb , d t Dj b = djDtb, 

where j = x, y, z. Denoting by u a the 20-component vector (Ab, Dtb, Djb), this system has the form of Eq. (JTJ) where 
the principal symbol 

A* emu? = (0,n j D jb , nj D tb ) 

is symmetrized by the trivial symmetrizer H a p = 5 a p. Therefore, the Cauchy problem for this system is well posed. 
If initial data is given for A b and its time derivative, both being divergence-free, then a unique solution exists, and it 
depends continuously on the given data. 

Next, let us consider the system on the half space E := {(x,y,z) £ R 3 : x > 0}. The normal matrix is obtained 
from the principal symbol with n = —dx the unit normal one-form. Therefore, the zero, in- and outgoing fields are 

D tb + D xb = (d t + d x )A b (incoming field) 

D tb - D xb = {d t - d x )A b (outgoing field) „ . _ „ . „ . „ . 

D yb = d y A b (zero speed field) dtAt ~ ° xAx + dyAy + dzAz ' 

D z b = d z Ab (zero speed field) 

where on the right-hand side we have explicitly written down the constraint equation in order to explain the problem: 
namely, when written in terms of the cigenfields, we sec that the constraint not only depends on the in- and outgoing 



1 See, for instance, Refs. [T3 | and [To| and references therein. 

2 Notice that no metric is needed to determine the one-form n c . 
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fields, but also on the zero speed fields! Therefore, imposing the constraint on the boundary in order to guarantee 
constraint preservation leads to a boundary condition which couples the incoming fields to outgoing and zero speed 
fields. 3 

Of course, this problem does not arise when writing Maxwell's equations in terms of electric and magnetic fields, 
which directly leads to a symmetric hyperbolic system with trivial constraint propagation. So the fact that the 
system is characteristic and has constraints with nontrivial propagation has conspired to convert a simple problem 
into an intractable one. Notice also that in the potential formulation there is no clear correspondence between the real 
radiation degrees of freedom and the in- and outgoing fields. Some of the fields describe gauge degrees of freedom. 

Until recently, the above problem was open. 

B. Well posed and strongly well-posed boundary conditions 

If one writes the equations as a first-order system, then at an initial surface the value of all the fields needs to be 
prescribed (in such a way that they satisfy the constraints if such are present) in order to obtain a unique solution. 
At the boundary the situation is different in general. Only data for the incoming fields must be given while the values 
the other fields acquire depend on the evolution in the bulk. Then, the question is whether or not one can control 
the fields in terms of the given initial and boundary data when imposing a boundary condition which couples the 
incoming fields with the outgoing and zero speed fields. Here, by "control" we mean the following: 

Definition 1 Consider a first-order strongly hyperbolic system of the form {I]) with initial data u Q | S[) = f a and 
boundary conditions of the form L r a (t,x,u)u a \^ = g r . This problem is called strongly stablefTH/ if there are norms 
II ' lis,. II ' II T bounding the fields and their derivatives on and T, respectively, and a constant e > such that each 
solution u satisfies an estimate of the form 

\\u\\ ST +s\\u\\r<F(T,\\fh ,\\g\\ r ), (2) 

where F(t,x,y) is a smooth function in its arguments that vanishes when t = or when (x,y) = (0,0). 

Remark 1 In particular, the estimate (0) implies the continuous dependence of the solutions on their data. A problem 
which is strongly stable and for which a unique solution exists on a suitable time interval [0, T] is called strongly well 
posed. 

A trivial example of a strongly stable problem is a symmetric hyperbolic system for which all the fields are incoming 
at the boundary, since in this case the values for all the fields must be specified there. More generally, a symmetric 
hyperbolic system with non- characteristic boundary surface and maximal dissipative boundary conditions @, | is 
strongly stable. A general theory for strictly hyperbolic first-order systems which allows to determine which boundary 
conditions lead to a strongly stable problem by purely algebraic methods was developed by Kreissfl7j. 

If a problem satisfies the same properties as in Def. [1] with the difference that the constant e in © is zero, then 
we say that the problem is stable. In this case, the estimate @ can still be used to show existence and uniqueness 
properties for some systems, but the set of systems for which this can be done is rather small, usually involving only 
linear systems or very restrictive boundary conditions. For this reason, in the following, we restrict ourselves to a 
smaller class of strongly hyperbolic systems for which strong well posedness can be shown. 

C. Second order systems 

A particularly interesting set of systems are those which are second-order, in particular those who describe systems 
of wave equations of the form 

g ab (<S> B )V a V b § A + F A (1> B , V c $ c ) = 0, A=1...N, (3) 



Instead of imposing the constraint itself on the boundary one might try to set some linear combination of its normal and time derivatives 
to zero, obtaining a constraint-preserving boundary condition that does not involve zero speed fields. Unfortunately, this trick only 
seems to work for reflecting boundaries, see Refs. Il5| and Il6(l for the case of general relativity. In our example, such boundary conditions 
are given by dtAt = d x A x = dtA y = dtA z = which imply dt(V b Ai > ) = 0. 
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where <& A is some set of tensorial fields, g ab a Lorentzian metric which in general depends locally on the fields, and V 
is a given connection on M . We impose initial conditions, 

* A k = *o> d t $X o =n$ (4) 

and boundary condition of the following form: 

N N 

(T d + aN d )V d $ A = Y, cAd B(t, x, $)V d $ B + Y dA s{t, *, ®)® B + G A (t, x), (5) 

B=l B=l 

where here a > 0, T is a future-directed time-like vector field tangent to T and N is the outward unit normal to 
T, see Fig. [TJ Furthermore, the matrix coefficients c Ad B have the property that they can be made arbitrarily small 
by an appropriate local linear transformation & A — > J A b& b of the fields (see Ref. [9| for the details). In particular, 
this is the case if the matrix operator c Ad b^J d is in upper triangular form with zeroes on the diagonal since then the 
matrix elements c Ad b can be made arbitrarily small by a simple rescaling of the fields, much like in the proof of the 
Lyapunov stability theorem. We have the following theorem: [j, H, @ 

Theorem 1 The IB VP (L2(/£J5p * s strongly well posed. 

This theorem constitutes the key result which allows the consideration of a whole new class of boundary conditions 
which are flexible enough to solve the IBVP of electromagnetism in the potential formulation and also the full 
Einstein vacuum field equations in harmonic coordinates. [3, This new result was obtained by realizing two things: 
first, that the above equations, when viewed as a pseudodiffcrcntial system 4 and then reduced to first-order is a non- 
characteristic system (see for instance Refs. [l8| and and second, that the Kreiss theory 17 applies equally well 



to pseudodifferential equations. The above boundary conditions are just of the type required in this theory to give a 
strongly well posed system. Later, the problem has been understood from a different point of view using more mundane 
energy estimates based on integration by parts instead of pseudodifferential calculus. @ The key observation is that 
for second-order systems like Eq. ([3]) there are many different hyperbolizations, and so many different energy norms. 
This freedom comes about because for a general Lorentzian metric there is no preferred time direction. Different 
time directions give rise to first-order reductions which differ among each others by constraint terms, and so they are 
incquivalent. By choosing the time direction in a suitable, but nonstandard form one obtains a non-characteristic 
first-order system which is amenable with the usual theory Notice that this situation is rather different than the one 
in fluid dynamics, where the four-velocity singles out a preferred time direction at each point of spacetime. 

D. Application: Maxwell's equations in the potential version 

As a first application, we return to Maxwell's equations for the four-vector potential on the half space E, 

V a V a A b = 0, V b A b = 0. 

The initial data A b \^ o = / b (0) , <9tA| Eo = / b (1) is subject to the conditions V h / 6 (0) = and V b / fc (1) = 0, so that the 
constraints are correctly propagated into the domain of dependence of the initial surface T,q. The boundary conditions 
are described in a geometrically elegant way by introducing the complex null basis of vector fields 

K:=d t ~d x , L:=d t + d x , Q:=d y + id z , Q:=d y -id z , i = \/-L 

Then, we consider the following set of conditions at the boundary T = [0,T] x <9£, 

^kA k = q K , (gauge) 

V k Aq = q' Q (Sommerfeld) (6) 

V K A L = -V L A K + VqAq + VqAq, (constraint V b A b = 0) 

where Ak '■= K b A b , := A" a V a , etc... and qK, are boundary data. The first two equations (which constitute 
three real conditions since Q is complex) are non-incoming conditions, usually called Sommerfeld conditions, since 



4 Think of this as the Fourier transform of the system. 
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they control the incoming fields. The left-hand side of the third equation is also an incoming field, but now it is 
sourced by first-order derivatives of the fields and Aq, for which the first two equations already give conditions 
to. Therefore, this system conforms with the hypothesis of Theorem [1] with a = 1, T = dt and N = —d x such that 
T + aN = K and leads to a strongly well posed problem. A simple way of understanding why this works for the 
second-order problem whereas it failed for the first-order reduction discussed in subsection III Al is the following: the 
first two equations in the boundary conditions ([6]) give rise to a strongly stable system for Ak and Aq, respectively, 
and hence Ak and Aq each satisfy an estimate of the form In particular, one controls all first derivatives of the 
fields at the boundary. Therefore, when applying the estimate to the wave equation for A^ with the third boundary 
conditions, the source terms on the right-hand side are controlled. 

Regarding the physical interpretation of the first two conditions in Eq. (O, the first controls the incoming gauge 
fields of the form A b = W b (e tL ' x ) while the second controls incoming electromagnetic fields of the form A b = Q b e lLx . 
However, a problem arises because V kAq is not gauge-invariant. For this reason, we replace the boundary condi- 
tions ([6]) by 

^ kA k = qx, (gauge) 

VkAq - V Q A K = q Q (radiation) (7) 

V K A L = -V L A K + VqAq + VqAq, (constraint V b A b = 0) 

where now the boundary data qQ represents the .fTQ-components of the electromagnetic field tensor which is gauge- 
invariant. The new boundary conditions (|7|) still satisfy the hypothesis of Theorem [T] and we obtain a strongly well 
posed problem. As we shall see, these new boundary conditions lead to gauge uniqueness, whereas the conditions (|6]) 
do not . 



E. Application: The linearized Einstein equations in harmonic coordinates 

As a second application we describe a strongly well posed IBVP for the Einstein field equations in vacuum when 
linearized on the flat background (M, -q) where M = [0, 00) X E, E denoting half space and rj — —dt 2 + dx 2 + dy 2 + dz 2 
the Minkowksi metric. The restrictions of linearization and working on half space are done for the sake of simplicity. 
A (local in time) well posed formulation for the full vacuum equations is given in Ref. H . When harmonic coordinates 
are used the linearized Einstein equations are 

V c V c h ab = 0, H a := V b (h ab - Uabri cd h cd ) = 0, 



where h ab denotes the first variation of the metric fields. As in the previous application, the initial data /i a f)| s = /lt!\ 
dth ab \-£ = f^l is chosen such that the constraints H a = and &tH a = are satisfied at Eo- In order to describe the 
boundary conditions, we use the complex null basis of vector fields K , L, Q, Q defined in the previous application. Let 
us start with the constraint conditions and work our way up on the triangular structure of the boundary conditions. 
The imposition of the constraints H a = at the boundary surface T yields 

^Kh Q Q = -V L,hKK + Vq/i K q + Vgh K Q, 

^Kh L Q = -^hh K Q + Vq/i^l + Vq/iqq, (8) 
Vjf/iiz, = -VLhQQ + V Q h L Q + VQh L Q. 

The right-hand side of the last equation contains first-order derivatives of the fields Iiqq and h^Q which are controlled 
in the first two equations. The right-hand sides of the first two equations require control of the first-order derivatives 
of hxK, h-KQ, hj<L and Uqq. Therefore, we impose the following gauge boundary conditions: 

^nh K Q - \^q^kk = qQ, (9) 
Vj<r/lKri - ' L h KK = q L , 

which control the fields hxK, hnQ and Kkl- Here, qn, qQ and qL are boundary data, and the purpose of introducing 
the second term on the left-hand sides of the last two equations will become clear in the next section. It remains to 
control the field /iqq, which is related to the field h ab = Q a Q b e %L x describing incoming gravitational radiation. Two 
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possibilities have been contemplated for controlling this field. The first specifies the shear, a := Q a Q b S7 a Kb, of the 
null congruence corresponding to the outgoing null vector field JsT.fl^-llII In the linearized setting considered here this 
yields the condition 

V K h QQ - 2V Q h KQ = 2d. (10) 

The boundary conditions (|9I10I8|) satisfy the hypothesis of Theorem [1] and therefore, one obtains a strongly well 
posed IBVP. 

However, in the next section we will see that the shear condition (|10p leads to difficulties when discussing geometric 
uniqueness. For this reason, it is convenient to replace this condition by a new condition which specifies the Weyl 
scalar ipo := —2K a Q b K c Q d R a b c d instead, where here, R a bcd denotes the linearized curvature tensor associated to h a b. 

This gives m m 

VI-Ziqq + V Q (V Q h KK - 2V K h KQ ) = Vo- (11) 

This condition involves second-order derivatives of the fields and handling it requires a little bit more work. However, 
using the pseudodifferential first-order reduction of Ref. [|| it is possible to prove that the boundary conditions (|9I11I8[) 
lead to a strongly well posed IBVP, sec Ref. [l9| . Estimates on the amount of spurious reflections introduced by the 
boundary conditions (fTO)) and (fTTj) were given in Refs. [HI, [24| and fl9| . 

III. GEOMETRIC UNIQUENESS 

In the previous section we have reviewed well posed IBVP for Maxwell's equations and the linearized Einstein's 
equations. These problems allow, from the point of view of partial differential equations, to construct unique solutions 
given appropriate initial and boundary data. In this section we discuss the question of gauge uniqueness, i.e. the 
relation that initial and boundary data need to satisfy such that the solutions they give rise to are gauge related. 

As discussed in the introduction, two initial data sets (h, k) and (h, k) on So are related to each other by a 
diffcomorphism of So if and only if their corresponding Cauchy developments on D(So) are related to each other by a 
diffeomorphism of £>(So) which leaves So invariant. It would be nice to have a similar statement for the IBVP, that is 
a similar statement that does not only consider two solutions on the domain of dependence D (So) of the initial surface 
but on the whole manifold M = [0, T] x S. Given two initial data sets (h, k) and (h, k) on So and two boundary data 
sets q and q on T := [0,T] x <9S satisfying suitable compatibility conditions at the edge S := {0} x <9S, one would 
like to know under which circumstances the corresponding solutions (M, g) and (M, g) are related to each other by 
a diffcomorphism on M which leaves So and T invariant. Is it, say, possible to relate q and q by a transformation 
on T alone, such that for given initial data on S the resulting metrics g and g are related to each other by a 
diffcomorphism? 

As pointed out in Ref. [25| there are several difficulties with this question: 

(i) It is a priori not clear what the boundary data q should be. Unlike for the case of the initial surface, q cannot 
represent the first and second fundamental forms of T, as explained in the introduction. 

(ii) The boundary data q = (qK,qQ,qL,&) and q = [qK , 9q, <?l, VOi respectively, introduced in subsection III El 
depend on a specific choice of the outgoing null vector field K, or at least on its orthonormal projection onto 
the boundary surface which is a time-like vector field T. It is not clear if in general there is such a preferred 
time direction on T. 5 

(iii) Related to the first two items is the question whether or not it is at all possible to find a transformation on T 
which relates q and q, independent of the initial data. 

One possibility for dealing with these issues is to introduce a background metric and to formulate the boundary 
conditions in a covariant way based on the covariant derivative defined by the background metric. This is the 
approach taken in Refs. [lj| and @ and further developed in Refs. [20| and [2l|. However, the issues above are then 
traded by the question on the dependency of the solution on the background metric. 

In the following, we analyze these issues in the much simpler setting of Maxwell's equations in the potential 
formulation, where the role of the diffcomorphism is played by the gauge transformations, and linearized gravity. 



5 However, see the recent proposal in Ref. |25| . 
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A. Gauge uniqueness for Maxwell's equations for the potential formulation 

We reconsider the source free Maxwell equation for the four- vector potential Ab on the half-space £, 

V a (V a A b - V b A a ) = 0. (12) 

In contrast to the previous examples, however, we do not necessarily impose the Lorentz gauge. The initial data 
for this problem are the three- vector potential and the electric field, {A^\e^), and the boundary data qQ which 
corresponds to the K Q-components of the electromagnetic field tensor and is related to the incoming electromagnetic 
radiation in the direction of L = dt + d x , see subsection III Dl Wc assume that this data is smooth, satisfies the Gauss 
constraint d 3 'Ej = and suitable compatibility conditions at S. Since the initial value of the magnetic field only 
determines A^ up to the addition of a gradient, we are looking for a solution of Eq. (|12p such that 

AA^ = Af + d jX (0 \ StAj - djM* = Ef\ V K A Q - V Q A K \ T = q Q> (13) 

for some smooth function x^ 011 ^o- The next result shows that this data determines a smooth solution Ab of 
Maxwell's equations which is unique up to a gauge transformation. 

Theorem 2 The IB VP U2113\) possesses a smooth solution Ab which is unique up to a gauge transformation Ab = 
A b + VbX- 

Proof: We first show the existence of a solution in the Lorentz gauge, for which Eq. ([12")) reduces to the system 
of wave equations V a V Q A b = 0. The initial data, (Ab\^ ,d t Ab\-E ), is chosen such that Aj\ s = Aj , d t Aj| s = 

Ej ^ + djAo\ Sg and 9jAo| So = d J A^\ where the initial data for Aq is chosen smooth but otherwise arbitrarily. 
Then, we evolve the well-posed IBVP with boundary conditions (|7J| described in subsection III Dt where the source 
function qx is chosen to be smooth and to satisfy the compatibility conditions at 5*, but otherwise is arbitrary. By 
construction, the resulting solution Ab solves V b Ab| Sn = 0, dtV b Ab\^ — 0, and therefore, it yields a solution of 

Maxwell's equations (|12|) in the Lorentz gauge with data (|13|) for x^ = 0. 

As for uniqueness, let A^ and A^' be two smooth solutions of Eqs. (|12I13|) . We show there exists a function 
X on M such that 8Ab := A^ — A^' = VbXi t na t is, the two solutions are gauge related on M. In order to 
see this, we first notice that SAb satisfies Maxwell's equations ([12)1 with initial data satisfying <5A/| E = djSx^ ' 
and dtSAj — dt5Ao\^ g = and trivial boundary data, qQ — 0. Next, wc transform SAb into the Lorentz gauge by 
finding \ sucn th & t V & V&x = ^ b SAb and setting Ab := SAb — V&%. The transformed potential Ab satisfies the 
Lorentz gauge. Furthermore, choosing initial and boundary data for x such that xls = $X : ^*xIe = ^-oIe > ano - 



Va-VxxIt = ^kA k \ t we obtain A t 



0, d t A t 



and VkAx 



= 0. Therefore, Ab is a smooth solution 

r 

of the IBVP described in subsection III Dl which has trivial initial data and satisfies the boundary conditions ()7J) with 
trivial data. By uniqueness of this problem it follows that Ab = 0, which implies that SAb = VfcX- rj 

Remark 2 From the proof of the theorem it follows that the choice for ^4ols an d the boundary data qx in the IBVP 
described in subsection UI D\ with boundary conditions corresponds to a pure gauge freedom. 

In contrast to this, the same IBVP with boundary conditions Q) replaced by the gauge- dependent condition 
\7kAq\j- = q'g, does not lead to gauge uniqueness. The reason is that a gauge transformation induces the transfor- 
mations qx = Qk + VkV^x an d q'g = Qq + Vx Vqx on the boundary data, which overdetermines the gauge function 
X- 



B. Geometric uniqueness for the linearized Einstein equations 

Einstein's field equations in vacuum, when linearized about the Minkowski metric rj a b, are 

- V c V c /i ab - VaVb(v cd h cd ) + 2V c V (a h b)c = 0, (14) 

where h a b denotes the first variation of the metric. The initial and boundary data are the first and second fundamental 
form of the initial surface, (fry , fcy '), and the Weyl scalar ip Q defined in Eq. (TTT)) . respectively. We assume the data to 
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be smooth, to satisfy suitable compatibility conditions at S, and to satisfy the linearized Hamiltonian and momentum 
constraints G ijrs didjh^ = and G^djkif = 0, where G ijrs := S^ r S s ^ j - 8^S rs . We arc looking for a solution of 
Eq. (Hi]) satisfying 

M So = hf) + 2d {l X n , d t h l3 - 2%^. )0 | So = -2(fcW + dftf), 

V 2 K h QQ + V Q {V Q h KK - 2V K h KQ )\ r = ifo, (15) 

where Xj and / arc a smooth vector field and a smooth function on Eo, respectively, representing the initial gauge 
freedom. 

Theorem 3 The IB VP \14]15)) possesses a smooth solution h a b which is unique up to an infinitesimal coordinate 
transformation h a b = h a b + 2V( a £j,) generated by a vector field £ . 

Proof: The proof works exactly as in the electromagnetic case, and is based on the harmonic IBVP described in 
subsection III El with boundary conditions (fTTj) . rj 

Remark 3 With respect to an infinitesimal coordinate transformation the boundary data qx, qQ and q^ introduced 
in subsection \II E] transforms according to 

q~K - qi< = 2V|-£if , q Q - q Q = V k ^q, q~L - qi = V^£l- 

The right-hand sides of these equations provide a complete set of boundary data for the gauge source vector £ a . 
Therefore, the choice for the data /tools > ^ojls an d Ik, 1q and qi, which is left unspecified by the physical data 

(h\j,k^ ) and ?/>o in* the harmonic IBVP corresponds to infinitesimal coordinate transformations. 

On the other hand, the shear a transforms as a = a — Vq£k ■ Therefore, replacing the boundary condition ill]) with 
the shear boundary condition tl(J\) over determines £ a at the boundary, and one does not obtain geometric uniqueness. 

Returning to the issues (i)-(iii) described above, we can say that in the simple case of linearization about Minkowski 
space we have solved the points (i) and (iii). That is, we have identified the boundary data for this problem that 
leads to geometric uniqueness. Regarding point (ii), the time-like vectors T that we have chosen in our problem 
correspond to the future-directed unit normal to the initial surface Eo at S, parallel transported along the geodesies 
orthogonal to Eo- It is clear that when considering the nonlinear case, or even for the case of linearization about a 
nonflat spacetime, several difficulties appear. These additional difficulties are mentioned in the next section. 

IV. CONCLUSIONS AND OPEN ISSUES 

In spite of our better understanding of the IBVP for Einstein's field equations from the point of view of partial 
differential equations, several problems are still open and they have mostly a geometrical character. Geometric 
uniqueness presents at least two main obstacles. One is the question of how to generalize the local-linearized version 
we have discussed above to the global, fully nonlinear case. The problem is that we have imposed no restrictions on 
the normal component £jv of the vector field generating the infinitesimal coordinate transformation. However, such a 
restriction is necessary in order to keep the boundary surface fixed under a diffcomorphism. Unfortunately, it does not 
seem possible to restrict £at with our current boundary conditions. As a consequence, the evolution of the boundary 
(viewed as an embedded surface in the unknown spacetime (M, g)) might depend on the initial gauge, say. Actually, a 
similar question arises when two solutions are considered on two different Cauchy surfaces and one asks the question 
whether or not the two solutions are related by a diffeomorphism. In order to answer this question, one needs to find 
the corresponding Cauchy developments. Our problem has this sort of difficulty also, and one could imagine solving 
this aspect of the problem " along the way" , that is while one is computing the time evolution, one can ask, step by 
step, whether or not a diffcomorphism relating the two solutions exist. This is what we call the zipper problem. One 
is mending the gap between the two boundaries as one is solving the evolution equations, see Fig. [2] 

The second obstacle is the choice of a time-like vector field T in our present boundary conditions. Contrary to our 
initial expectation this problem might not be of a fundamental character; it is rather a manifestation of our inability 
to specify a non-incoming radiation condition correctly. Indeed, our present conditions single out the outgoing waves 
with wave vector along the null direction K we have chosen. In this sense, our boundary condition is perfectly 
absorbing for such outgoing waves, but not for outgoing waves with different wave vectors. A genuine non-incoming 
wave boundary condition should be independent of any specific null or time-like direction at the boundary, and can 
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only depend on the normal vector to it. For simpler systems like the scalar wave equation this is indeed the case: [26| 
the imposition of perfectly absorbing boundary conditions leads to a nonlocal condition which is independent of a 
preferred time direction at the boundary. Although this condition might not be practical because of its nonlocality, it 
is possible to derive a hierarchy of local boundary conditions depending on a null vector K approximating the nonlocal 
condition with better and better accuracy. So it is expected that the choice of K becomes less and less relevant as 
one moves down the hierarchy. It should be interesting to explore this idea for the case of general relativity, departing 
from a nonlocal boundary condition which is independent of K . 6 

We end this article by asking whether or not we can apply the methods discussed here to metric formulations of 
Einstein's field equations other than the harmonic one. For example, a formulation that is widely used in numerical 
calculations is the BSSN formulation [H, [29[ which yields a mixed first/second order strongly hyperbolic evolution 
system]?! [13, HU- There are two ways of applying the theory discussed in this article to this problem. The first 
consists in finding a pscudodiffercntial first-order reduction of the evolution system which is non-characteristic. Such 
a reduction should be enough to impose constraint-preserving boundary conditions which are strongly stable, and for 
which the Krciss theory can be applied. The second way is to exploit the non-uniqueness of time directions in general 
relativity and to write down a symmetric hyperbolic first-order reduction with respect to a suitable time evolution 
vector field, such that the boundary surface is non-characteristic and constraint-preserving boundary conditions can 
be specified in maximal dissipative form. We think that both ways should, in the end, be equivalent but this needs 
further investigation. Regarding the BSSN system it turns out that it possesses zero speed fields which are intrinsic to 
the formulationQ, and therefore, it does not seem possible to apply the above methods unless one considers moving 
boundaries. Nevertheless, partial results[32[ have been recently obtained based on different methods. 
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